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Summary 

Numerical comparisons are made of an exact 
solution with the Eikonal approximation to the 
Lippman-Schwinger equation for an equivalent first- 
order optical potential for heavy ion scattering. The 
exact solution proceeds by using a partial-wave ex- 
pansion of the Lippmann-Schwinger equation in mo- 
mentum space. Calculations are made for the to- 
tal and absorption cross sections for nucleus-nucleus 
scattering. The results are compared with solutions 
in the Eikonal approximation for the equivalent po- 
tential and with experimental data in the energy 
range from 25 A to 1000 A MeV. The percentage dif- 
ferences between the exact and Eikonal solutions are 
shown to be small above 200 A MeV, but at lower 
energies, differences become large and corrections to 
the first-order optical model are required. 

Introduction 

In high-energy nucleus-nucleus scattering, many 
nucleons can interact mutually and the structure 
of multiple scattering is richer than nucleon-nucleus 
scattering. A simple approach in nucleus-nucleus 
scattering is to consider the scattering in terms of 
each constituent of the projectile nucleus interacting 
with each constituent of the target nucleus. Other 
terms may contribute to the scattering, such as a pro- 
jectile constituent interacting consecutively with two 
different constituents of target (i.e., double scatter- 
ing). Similarly, contributions may come from three, 
four, and more successive scatterings. The formal- 
ism using this approach is called multiple scattering 
theory. 

The nucleus-nucleus scattering processes are con- 
veniently analyzed by employing an optical potential 
theory. Once the optical potential is determined, the 
original many-body scattering problem reduces to a 
two-body scattering problem. However, the price of 
reducing a many-body problem to a two-body situa- 
tion is that the optical potential will be, in general, 
a complicated nonlocal, complex operator. Thus, for 
practical applications we shall require the approxi- 
mation method to determine the optical potential. 

An early generalization of the optical model ideas 
in nuclear physics to the study of alpha-decay of 
nuclei was made by Ostrofsky, Breit, and Johnson 
(ref. 1). Bethe (ref. 2) introduced the concept of an 
optical potential in order to describe low-energy nu- 
clear reactions within the compound nucleus model. 
The description of high-energv nuclear collisions by 
means of the optical model formalism was initiated 


by Fernbach, Serber, and Taylor (ref. 3), who first 
tried to describe elastic nucleon- nucleus scattering 
in terms of nucleon- nucleon collisions. They argued 
that at high energies, a nuclear collision should pro- 
ceed by way of collisions with individual target nucle- 
ons by using known nucleon-nucleon cross sections. 
This multiple scattering analysis led to the conclu- 
sion that particles should move more or less freely 
through nuclear matter at high energies. The fact 
that the optical potential is complex is worth notic- 
ing. The imaginary part of the optical potential cor- 
responds to the absorption of the incident beam by 
target nuclei, and the real part of the optical poten- 
tial corresponds to the refraction of the beam with- 
out any disturbance to the target nuclei. Watson 
(ref. 4) and Kerman, McManus, and Thaler (KMT) 
(ref. 5) developed the formal theory of the scatter- 
ing of high-energy nucleons by nuclei in terms of the 
nucleon- nucleon scattering amplitude. 

A general multiple scattering theory for the scat- 
tering of two composite nuclei (neglecting three-body 
interactions) has been developed by Wilson (refs. 6 
and 7). Calculations of the reaction and absorption 
for the heavy ion projectile was well developed by 
Wilson and Townsend (refs. 8 and 9) by using an 
Eikonal approximation to a first-order optical model. 
In the first-order optical model, the excitation of the 
projectile or target in intermediate states is neglected 
in elastic scattering. The Eikonal approximation is 
based on a forward scattering assumption, as well as 
considerations of the strength of the potential (refs. 6 
and 10). A second-order solution to the Eikonal 
coupled-channels (ECC) model was developed by 
Cucinotta et al. (refs. 11 and 12) and was found to 
give improved accuracy over the first-order solutions 
in limited studies for several collision pairs and en- 
ergies. The Eikonal approximation is computation- 
ally efficient because it requires only a few numeri- 
cal integrations for implementation. In this report, 
neutron- nucleus, alpha- nucleus, and carbon-nucleus 
total and absorption cross sections are calculated by 
using the Lippmann-Schwinger equation. By com- 
paring calculations with identical physical inputs 
that use the exact solution of the first-order opti- 
cal model to the Eikonal approximate solutions, we 
provide an important validation of databases used in 
cosmic-ray jstudies. The Eikonal model is also un- 
able to account for nuclear medium corrections, al- 
though such studies may be performed in the future 
by using momentum space methods. The treatment 
of medium corrections will be required if further im- 
provements in nuclear databases are needed. 



In the remainder of this paper, we present 
the formalism for the first-order optical model and 
partial-wave decomposition of the Lippman- 
Schwinger equation. In the following sections, we dis- 


cuss the optical potential, the Eikonal model, model 
calculations describing physical inputs, and the com- 
parison of exact and Eikonal solutions with experi- 
mental data. 


First-Order Optical Model and Partial- Wave Decomposition of Lippmann-Schwinger 
Equation 

The total Hamiltonian H satisfies the time-independent Schrodinger equation 

W) = E\i>) (1) 

where E is the eigenvalue of operator H and |t/>) is the state vector. In case of nucleus-nucleus scattering, the 
total Hamiltonian is given by 

H = H P + H t + V (2) 

where V is the potential energy operator between projectile (P) and target (T) constituents. Here, V is defined 
as 

A P A t 

v =a o) 

j = i q=i 

where j denotes the projectile, a denotes the target nucleus, v aj is the interaction potential, and Ap and Ap 
denote the projectile and target mass numbers, respectively. 

The transition operator for scattering the a constituent of the target with the j constituent of the projectile 
is defined by the Lippmann-Schwinger equation (which is the integral form of the Schrodinger equation) as 

T a] = V aj + V aj G 0 T Qj (4) 

(ref. 6) where G 0 is the Green’s function given by 

G 0 = (E-Hp - Hr)- 1 ( 5 ) 

The solution of the multiple scattering problem is generally intractable in the exact form. An approximate 
potential is developed by using the transition operator of equation (4) to replace the highly singular two-body 
potential. This is usually developed through an optical-potential formalism following the methods of Watson 
(ref. 4) or Kerman, McManus, and Thaler (ref. 5). Wilson and Townsend (refs. 6-8) have considered a first-order 
optical potential in the impulse approximation for nucleus-nucleus scattering which neglects the excitation of 
the projectile and target in the intermediate states for elastic scattering (coherent approximation). Here, the 
transition operator (T) obeys 

T = U + Ug 0 T (6) 

where g Q is the free two-body Green’s function and the optical potential is given by 

u = £ t Q} ( 7 ) 

a 7 

where 

tatj — u otj + u aj9o^aj ( 8 ) 
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An alternative derivation of a first-order optical potential for nucleus-nucleus scattering that considers 
antisymmetrization effects following the approach of reference 5 would lead to 


u = z 1 t 

ApA T 


(9) 


which is approximately the same as equation (7) for ApA T » 1. 

To solve equation (6) we take the matrix element with respect to the projectile and target ground states. 
If k j and k 2 are wave vectors of two nucleons in the lab system, we define the relative wave vector as 


k = 


ki - k 2 


If we take the half on-shell matrix element of T in equation (6), we get 

<k'|T|k) = <k'|U|k) + f ( k 1 u l k ")( k "l T l k )< fl{ " 

j E k - E k n + IT] 

where |k) is the on-shell wave vector. The R matrix satisfies 

(k'lRjk) = (viuik) + p [ vvwrmM' 

J E h - Ey> 


( 10 ) 


( 11 ) 


( 12 ) 


where P denotes the principal value of the integral. The R matrix is related to the T matrix by the Heitler 
equation, which is given by 

T = R-urR,5(£ , -/i 0 )T (13) 

after the partial-w r ave decomposition 

Rlik.k') = U t (k, ~ ~ ~ k " dk " (14) 

o 

where p is the reduced mass. Thus, we have used the following expression for the partial wave decomposition: 

(15) 


( k/ |R|k) = J2 k) P]{k\ k) 


The on-shell 7} is related to the on-shell R] via the Heitler equation. Thus, 

rp _ Rl(k, k) 


1 T A£kR](k : k) 


( 16 ) 


The nucleus-nucleus transition matrix T(g) and the scattering amplitude f(q) for the nucleus-nucleus system 
are given by 


T (g) = T i p '( cos 


1=0 


and 




(17) 


( 18 ) 
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respectively, and the differential cross section is given by 

% - <i9) 

Now, by integrating over a solid angle we obtain the total elastic cross section as 

<7el = / \f(q)\ 2 dQ (20) 

The total cross section is related to the forward scattering amplitude by the optical theorem as 

tftot = X Im /el(0 = 0) (21) 

where Im denotes the imaginary part. Thus, after T) is calculated from equations (14) and (16) by the matrix 
inversion method, the scattering amplitude can be obtained from equation (18) and the total cross section can 
be obtained from equation (21) by taking the imaginary part of the forward scattering amplitude. Now, the 
absorption cross section is simply obtained by 


G r — Otot &e\ 


( 22 ) 


The expression for the scattering amplitude given by equation (18) requires an infinite number of partial 
waves for the t matrix to be calculated from equation (16), and usually truncating the sum to a finite number 
of partial waves is necessary. The results obtained are very reasonable at low energy. It is well known that 
as the projectile energy increases, the number of partial waves necessary also increases to make the sum 
in equation (17) divergent. It is also well known that for high partial waves and high energies, the Born 
approximation for the 1th partial-wave component of T/ is a good choice (i.e., T) as U{). We will utilize this fact 
to our advantage in order to make a correction for the contribution of higher partial waves. First, we rewrite 
equation (17) as 

^max o) i 1 00 11x1 

T(g) = £ =±i T/P/ (cos#) + £ -±- T; Pi (cos #) (23) 

/=0 ^max + 1 

where l m3X + 1 is the partial wave in which the Born approximation is good. This can be easily done by 
comparing C/j and the calculated 7} for each partial wave while solving equation (14). When making the Born 
approximation for the second term of equation (23), we obtain 


(max oi i i °° 9/11 

T( 9 )=£ T;P; (COS #) + Y, -±-UiPt(C0s6) 

1—0 


lu 


In 


i max 07 t 1 t max 9/ i 1 

+ E u ‘ P,icos "Ttt cos 0) 

1=0 1=0 


(24) 


in which we have added and subtracted a term in the last equation. Now, the second and third terms can 
combine to give us the three-dimensional U (<?), because it is a sum of all partial waves for U. Therefore, we 
finally get 

(max 97 1 -I (max 9/ 1 1 

T(g) = £ T,fl(. cos #)-£ — £ f/,P,(cos 6) + U( 9 ) (25) 

7=0 ^ 7=0 ^ 

Equation (25) is the final result, which tells us that in order to obtain the contributions from the higher partial 
wave, we have to calculate 7) only up to a certain Ixnax &t which the Born approximation becomes good. Then, 
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T(q) is obtained by summing 7} to l max and U (q) is added. At the end, the sum of U t up to l max has to be 
subtracted to avoid double counting. The phase shift 6/ is given bv 


T t = 


-h 2 e 2i6 ‘ - 1 


7rp 2 ik 

The total elastic and absorption cross sections are defined (ref. 10) as 


CT el = pE( 2/ + 1 )K e2 * Re *'- 1 | 2 


1=0 


and 


= j^E ( 2/ + !) (1 - rf) 


1=0 


respectively, where Re denotes the real part and 


*1 


_ e -2 Im 6 t 


Then, the total cross section is defined as 


(26) 


(27) 


(28) 


(29) 


2n 

atot = 


00 


E 


(21 + 1) [1 - rji cos (2 Re 6/)] 


(30) 


Optical Potential 

In the first-order optical model in the impulse approximation, the optical potential is the matrix element 
over target and projectile ground states of free two-body amplitude. Thus, 

</ |u N> = — A 4 p 7 A ~ 1 <0p0 r k'|£mo P o r k> (31) 

P T CCj 

At this point we digress to discuss and recall that the optical potential has a spin dependence that 
arises from the spin dependence of the t matrix. From symmetry principles, one can write the nonrelativistic 
t matrix for nucleon-nucleon scattering in terms of Pauli spin operators as 

t =J 4 + C (a i + cr 2 ) • h + M a\ ■ mer 2 • m + G a\ ■ h <r 2 . fi + Hcr\ l rr 2 • / 

+ D(a i ■ m (T2 ■ l + cry ■ l ct 2 ■ m) ( 32 ) 

(ref. 13) where cr is the Pauli spin operator and (/, m, n) defines a right-hand coordinate system. 

If we consider only spin-zero nuclei, and the spin projection of target and projectile nucleus is integrated 
out, the terms of t that are linear in the projectile and target nucleon spin are vanishing and leave only the 
first term A. 

Thus, we get the optical potential for spin zero-spin-zero, nucleus-nucleus (NN) scattering as 

ApA "T — 1 f 

(/|U|l) = ~ ApAx J dq Mq) Fp{q) ° T{q) ( 33 ) 

where F P is the projectile form factor and Gj is the target form factor. Therefore, we write 


A(q) = 


-h 2 

(2^)2p 


f NX (q) 


(34) 
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Eikonal Model for Total and Absorption Cross Sections 

The generalization of the Eikonal approximation to many-body scattering problems is given by Glauber 
(refs. 14 and 15), who applied it extensively to high-energy hadron-nucleus scattering. Wilson (ref. 6) 
has discussed the Eikonal approximation for the nucleus-nucleus optical model by using a coupled-channels 
formalism. In the Eikonal coupled channel (ECC) model (refs. 6 and 11), the matrix of scattering amplitudes 
for all possible projectile-target transitions is given by 

f(q) = j d 2 b e iq b {e l *< b > - l} (35) 

where barred quantities represent matrices, b is an impact parameter vector, q is the momentum transfer 
vector, and k is the projectile target relative wave number. In equation (35), Z is an ordering operator for 
the Z-coordinate, which is necessary only when noncommuting two-body interactions are considered (ref. 11). 
The matrix elements of x(b) are written as 

(0p0p|x(b)|0p0 r ) = ^ ^ Y : J d 2 qe l(i ' h F P {-q) Gp{q) / NN (g) (36) 

aj 

The first-order approximation to the elastic amplitude, which is obtained by neglecting all transitions between 
the ground and excited states and assuming that transitions between excited states are negligible, leads to 

fW( q ) = J d 2 b e- iq b {e^ - l} (37) 

(refs. 6 and 11) where 

X(b) = / d 2 q e i q b F P (-q) G T {q) f m (<?) (38) 

The total cross section found from the elastic amplitude by using the optical theorem from equations (21) 
and (37) is given as 

00 

o tot = 471- J bdb {1 - exp[(-Imx)cos(Rex)]} (39) 

0 

The total absorption cross section that is found by subtracting the total elastic scattering cross section from 
equation (39) is 

oo 

a r = 2Trj bdb {l— e- 2Im *} (40) 

0 


Model Calculations 

The two-body scattering amplitude is parameterized as 

/,mfa) = £ ^*NNe- s « ! '' 2 (41) 

where a is the nucleon-nucleon total cross section, a is the ratio of the real part to the imaginary part of the 
forward two- body amplitude, and B is the slope parameter. 

The one-body form factor is written in terms of the charge (CH) form factor as 


F{q) = 


■FcHfa) 

F p (q) 


(42) 
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where F p (q) is the proton form factor. For light nuclei (A < 16), we use the harmonic well distribution 


^CH = {l~sq 2 )e aql 


( 43 ) 


where the values for the parameters s and a are taken from reference 9. For medium and heavy mass nuclei 
(A > 17) where a Wood- Saxon density is appropriate, 


>) = 
^CH 


Po(r) 


l-exp(^) 


(44) 


We use an exact Fourier transform from equation (44) to obtain the charge form factor found in the series 
solution given by 

47 r 

F cn(q) = — po{q) 4>(q) (45) 


(ref. 16) where 

<j>(q) — nRc < 


— cos {Rq) 7TC sin (Rq) coth(7r cq) 
sinh(7rc^) R sinh(7rcg) 


2c , . m meg e mfl / c 
ttR [(c#) 2 + m 2 ] 2 


m= 1 


(46) 


The series in equation (46) converges rapidly, and the first three or four terms are accurate for most applications. 
Values for the parameters c and R are taken from reference 9. 


Discussion 

By using the formalism described in previous sec- 
tions, total and absorption cross sections for neu- 
tron, He, and C colliding with different target nuclei 
have been calculated in an energy range from 25 A 
to 1000A MeV. By using the Lippmann-Schwinger 
calculation and the Eikonal model, theoretical pre- 
dictions for total and absorption cross sections are 
compared in figures 1-6 with representative experi- 
mental data (refs. 18-22). All physical inputs (form 
factors and tw^o-body amplitudes) are kept identical 
in the Lippmann-Schwinger and Eikonal model cal- 
culations. The agreement is excellent between the 
Lippmann-Schwinger calculation, the Eikonal model, 
and the experimental data at higher energy for total 
and absorption cross sections. 

Calculations of the total cross sections for the 
nucleon- nucleus, helium-nucleus, and carbon-nucleus 
systems are shown in figures 1, 2, and 3, respec- 
tively. We observe from our calculations that at a 
lower energy, the percentage differences between the 
Eikonal model calculation and the exact calculation 
are higher when compared with those at higher en- 
ergy. This is an indication that the Eikonal model 
prediction for scattering cross sections is fairly accu- 
rate at higher energies. The Eikonal model calcula- 
tions are consistently lower than the exact calcula- 
tions because of the forward scattering assumption 
of the Eikonal approximation. Although the scatter- 
ing is dominated by forward angles at high energies, 


some contribution from large-angle scattering is al- 
ways present and is not represented by the Eikonal 
model. We also observe that both the Eikonal model 
and the exact calculations are well within the range 
of experimental data. At low energy, an improvement 
in the calculations will most likely occur by consid- 
ering the corrections to the impulse approximation, 
correlation effects (refs. 17 and 23), or perhaps rela- 
tivistic effects (ref. 24). 

In figures 4, 5, and 6 we present compar- 
isons of the partial-wave solution for the absorption 
cross section to the Eikonal model solution for the 
nucleon-nucleus, helium-nucleus, and carbon- nucleus 
systems, respectively. Experimental data are shown 
in these comparisons when available. The Eikonal 
model is seen to represent the exact solution quite 
well above projectile energies of 200A MeV, but be- 
low 200A MeV the differences are large. In some 
cases, the Eikonal model represents the experiment 
better than the exact solution, which is a definite 
indication that the first-order optical potential does 
not represent the complete physics at these lower en- 
ergies because the Eikonal model is an approximation 
to the exact solution for the equivalent potential. Ta- 
ble 1 presents a comparison of the experimental data 
from reference 20 with the Eikonal and Lippmann- 
Schwinger calculations for carbon- nucleus absorp- 
tion cross sections at an energy of 83.4 MeV. Al- 
though the calculations agree satisfactorily with ex- 
periment, further investigations of optical potential 


7 



theory, including nuclear medium effects, will be re- 
quired for the theoretical evaluation of absorption 
cross sections with high precision at lower energies. 

Table 1 . Comparison of Experimental Data a With 
Calculations for Carbon-Nucleus Absorption 
Cross Sections 


[■^lab / nucleon = 83 A MeV/ nucleon] 


System 

Values of a r , mb, for — 

Experimental 

Lippmann-Schwinger 

Eikonal 

12 C + 12 C 

960 ± 30 

874 

849 

12 C + 27 Al 

1400 ± 40 

1419 

1477 

12 C + 40 Ca 

1550 ± 60 

1750 

1737 

12 C + “Fe 

1810 ± 100 

2123 

1997 


a For the carbon nucleus, the experimental data are taken from 
Kox et al. (ref. 20). 


Figure 7 shows the number of partial waves re- 
quired to calculate the total, absorption, and elastic 
cross sections at energies of 100 A and 1000 A MeV 
for nucleon-carbon and carbon- carbon systems. We 
observe from the calculations that if we increase the 
energy of the projectile, we will need more and more 
partial waves for the calculation of cross sections. 
Figure 8 shows the total and absorption cross sections 
as the function of the slope parameter at energies of 
100A and 1000A MeV for the n- 16 0 system. We 
observe from our calculations that because the total 
and absorption cross sections saturate at large values 
of the slope parameter, a limiting value is reached. 

Concluding Remarks 

The exact calculation of the Lippmann-Schwinger 
equation for a first-order optical potential was con- 
sidered for light and heavy nuclei-induced reactions 
in an energy range from 25 A to 1000 A MeV. Com- 
parisons with Eikonal approximate solutions and ex- 
perimental data for target nuclei most important 
for space radiation studies were quite favorable and 
provided validation for earlier calculations used in 
cosmic-ray transport codes. The percentage differ- 
ences at high energy between the Eikonal calcula- 
tion and the exact calculation were small compared 
with those at low energy. Further improvements in 
the calculations will most likely occur by considering 
corrections to the impulse approximation and corre- 
lation effects, rather than to the first-order optical 
potential model employed in this paper. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
January 6, 1995 
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(a) n- 4 He scattering. 



(b) n- 12 C scattering. 


(e) n- 64 Cu scattering. 


(f) n- 208 Pb scattering. 


Figure 1. Comparison of total cross section calculation using Eikonal approximation with exact solution using 
equation (30) for nucleon- nucleus system in energy range from 25 A to 1000A MeV. Available experimental 
data are shown by error bars. 
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Figure 2. Comparison of total cross section calculation using Eikonal approximation with exact solution using 
equation (30) for helium-nucleus system in energy range from 25 A to 1000^4 MeV. 



Figure 3. Comparison of total cross section calculation using Eikonal approximation with exact solution using 
equation (30) for carbon-nucleus system in energy range from 25A to 1000^1 MeV. 
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(f) n- 208 Pb scattering. 


Figure 4. Comparison of absorption cross section calculation using Eikonal approximation with exact solution 
using equation (29) for nucleon-nucleus system in energy range from 25A to 1000A MeV. Available 
experimental data are shown by error bars. 
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Figure 5. Comparison of absorption cross section calculation using Eikonal approximation with exact solution 
using equation (29) for helium-nucleus system in energy range from 25 A to 1000 A MeV. Available 
experimental data are shown by error bars. 



Figure 6. Comparison of absorption cross section calculation using Eikonal approximation with exact solution 
using equation (29) for carbon- nucleus system in energy range from 25 A to 1000 A MeV. Available 
experimental data are shown by error bars. 
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(a) 100,4 and 1000A MeV for n- 12 C system. 



(b) 100A MeV for 12 C- 12 C system. 

Figure 7. Comparison of total, absorption, and elastic cross sections as a function of number of partial waves. 
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